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If QCD is to undergo a second order phase transition, the 
light quark return probability is universal for large times at 
qq ■ the critical point. We show that this behavior is distinct from 
' the one expected at the mobility edge of a metal-insulator 
, transition or a percolation transition in d< 4. Our results are 
accessible to current lattice QCD simulations. 

PACS numbers: 11.30.Rd, 12.38.Aw, 64.60.Fr 



In section 4 and 5, we contrast our results to the expected 
ones from a metal-insulator and a classical percolation 
transition. In section 6, we use semi-classical arguments 
to show what this means for the light quark spectrum. 
Our conclusions are in section 7. 



II. QUARK RETURN PROBABILITY 
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I. INTRODUCTION 

Light quarks trapped in a finite Euclidean volume V 
behave like electrons in small metallic grains The 
spontaneous breaking of chiral symmetry amounts to dif- 
fusing quarks in proper time with a vacuum diffusion con- 
stant D w 0.22 fm [0. For times larger than the ergodic 
time r crg = \fV / D or small virtualities, the quarks un- 
dergo ergodic motion, while for times shorter than the 
elastic time = l/2mg ~ 0.33 fm H] or large virtuali- 
ties their motion is ballistic. Here ttiq is the constituent 
quark mass in the vacuum. Diffusion sets in at interme- 
diate times and away from the quantum regime with the 
Heisenberg time tjj = 1/A = HV/ir and £ the quark 
condensate. 
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The concept of the quark return probability in the 
QCD vacuum as a chirally disordered medium (!]], bor- 
rows on the concept of the electron return probability 
in disordered metallic systems as first introduced by An- 
derson H in the context of localization. For long times 
and in the presence of sufficient disorder, Anderson noted 
that the electron return probability is finite. A similar 
observation can be made for the quark return probability 
in the QCD vacuum M. 

In this paper we would like to show that at the critical 
point of a second order phase transition, a quantitative 
change in the character of the disorder takes place, with 
consequences on the long time behavior of the quark re- 
turn probability. In section 2 and 3 we analyze the quark 
return probability in the vacuum and at the critical point. 



The eigenvalue equation of the massless Dirac operator 
for quarks in the fundamental representation and in the 
fixed gluon field A 



lf[A] q k = X k [A] q k 



(1) 



allows us to extend the theory into 4+1 dimensions with 
proper time t, and define the probability p(t) for a light 
quark to start at x(0) in V and return back to the same 
position x(t) after a duration t, as 

p(t) = ^(|(x(0)|e i ™+^)l t l|a ; (0))| 2 ) 4 . (2) 

The operator if[A] in (]l|) acts as a four-dimensional 
Hamiltonian for the evolution in proper time t. Indeed, 
the expectation value in (Q) is that of the proper time 
evolution operator, with the background gluon field A 
acting as a t-independent (static) random potential. The 
averaging in is over all gluon configurations using the 
unquenched QCD measure, with light quark flavors of 
current mass m (flavor symmetric case). The normaliza- 
tion in (|^) is per state, where N is the mean number of 
quark states in the four-volume V. 

In terms of the quark-eigenfunctions (jl]) in V, the re- 
turn probability reads 
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where the exponent e~ 2m l*l is solely due to the current 
quark mass in (^). We note that (||) is gauge- invariant 
and amenable to lattice QCD simulation. This form is 
best suited for numerical estimates. 

For analytical considerations, it is best to rewrite (||) 
in terms of the standard Euclidean propagators for the 
quark field, 
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p(t) = Slim f '1±!±L< • 



iV J (2tt) 2 

x (Tr(S{x,y; Zl )S\x, y,z 2 ))) A (4) 



with z\ 2 = m — iAi 2, and 



S(x,y;z) = (x\ 
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(5) 



Setting A12 = A ± A/2 and neglecting the effects of A in 
the averaging in (Q) (this is certainly true near the zero 
virtuality point), we find that in the flavor symmetric 
limit, the correlation function in (Q) relates to the pion 
correlation function for a proper analytical continuation 
of the current quark mass M . Specifically, 
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with z = m — iX/2 and E = J dA. For 2 = to, pion-pole 
dominance (long paths) yields 



Cn(x, y;m) « — ^ 
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with = rif/lir/L in V = L 4 and S = 

Using the GOR relation F 2 m 2 = mS, and the analyt- 
ical continuation m — > m — iA/2, we find [jjj 



C^ar, y;z) w — ^ 
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with the diffusion constant D = 2F 2 /T,. Inserting (||) 
into (|), and noting that E/A = N and p = 1/AV, 
where A is the mean quantum level spacing with £ = 
7T/9 (Banks-Casher relation), we conclude after a contour 
integration that 



p(t) = e 
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The validity of ( |10| ) is for t c i < t < tjj- For miff ~ 
ml/ <C 1, ( P-0|) is dominated by the constant pionic mode, 
giving p(t) ~ 1. What happens to the present arguments 
if QCD is to undergo a phase transition? 



III. SCALING AND UNIVERSALITY 

During a second order transition, the quark conden- 
sate and the pion parameters undergo structural changes 



which are governed by scaling arguments 
(||) takes the general form 
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where is the pion wavefunction renormalization, and 
Xn the pion susceptibility. At the critical point Tq 
(generic of a second order transition), V — L 3 /Tc- The 
representation (11) holds for long paths, or Euclidean 
momenta of order A < (2irTc). For larger Euclidean mo- 
menta or short paths, ( |ll| ) is saturated by quarks in the 
dimensionally reduced theory || . 
At the critical point, we have 



Xtv = S 
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where j3, v, S, rj are critical exponents. The prefactors in 
( |l2| ) are dimensionful constants that relate to the equa- 
tion of state (the last equation). For convenience they 
were set to 1. Hence, 
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where A = B - 1 + 1/6 and B = 4/(5(2 + 77) after using 
scaling relations. For a mean-field transition there is no 
wave-function renormalization and 77 = 0. Therefore: 8 = 
3, A = and B = 2/3. After the analytical continuation 
to — * to — iX/2, we obtain 



C (x vz)~ 1 V c «H*-») ( m -^/ 2 ) A 
^x,y,z)~ > e Q 2 + (m _ iA/2)B - 

As a result, the quark return probability becomes 
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(15) 



where C = VE/N2 A - B . 

In contrast to the vacuum case (^), the result ( |l4|) at 
the critical point displays new singularities. For fixed Q 2 , 
the pole in the A-plane in (^) characteristic of the vacuum 
phase is now changed to a branch point at A = — 2im, 
and a set of poles (Q 7^ 0) 



An 



-2im + 2ie 



i(2ri+l)ir/B|Q2| 1 /- B 



(16) 



with n = 0,1, ...,n m ax, and n max is the number of unit 
roots to 1 + z B = 0. The contribution of the cut is 
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(17) 
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while the contribution of the poles is 



IV. METAL-INSULATOR TRANSITION 
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The primed sum in ( J18| ) retains only those poles in the 
lower half of the A-plane. The result is still real as the 
poles occur in (A, —A*) pairs. 

For large times, the dominant contribution to p(t) re- 
sults from the Q — part (zero mode) of (fl7|). Setting 
the Heisenberg time to be t H = 1/A* - y^/U+'J) g ; 
and the mean number of levels to be N — E/A* ^> 1, 
the quark return probability simplifies 
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For a* = 2mV 5 /( 5+1 ) <C 1 (ergodic regime), the return 
probability for large times t is universal at the critical 
point with p(t) ~ {tn/t) 1 ^ 6 ■ This behavior is to be con- 
trasted with the vacuum result p(t) ~ 1 in_the ergodic 
regime and (T eie /t) 
We note that ( |14| ) 
tion in d=4 at m = (with the cutoff A) 



^ 2 in the diffusive regime [1J. 
obeys an 'anomalous' diffusion equa- 



-£>(A)V 2 + C v (x,y;-iX) = 



(20) 



for x ^ y. At the critical point the diffusion constant 
is A dependent and complex D(\) = 2 B (— iA) 1_B . For 
A ~ m ~ (zero virtuality) the diffusion constant van- 
ishes as |D(m)| ~ to 1_s . In this case, the Heisenberg 
time is tn = 1/A* ~ V s ^ s+1 \ and the ergodic time is 
r crg = 1/E C ~ /m}~ B , where £? c is the Thouless en- 
ergy Q . Since the constituent quark mass (half the sigma 
mass) becomes comparable to the pion mass at the crit- 
ical point, the elastic time is Td ~ \jm B l 2 (see (fl3|)). 
To assess the hierarchy of scales near the critical point, 
we enhance artificially the contribution from the ergodic 
regime by setting m/A* ~ 1 in power counting, in gen- 
eralization of a previous argument ||. For a mean- field 
transition: Td ~ V 1 ^ 4 and r crg ~ tH ~ F 3 / 4 . Hence the 
Ohmic conductance <tl = tH/T crg ~ 1 a situation remi- 
niscent of the metal-insulator transition, except that in 
the present case the density of states vanishes as A 1 / 5 . 
The time scales are still ordered in the thermodynamical 
limit, with the diffusive regime stretching all the way to 
the quantum regime. 

We note that the same scaling arguments we have used 
at T = T c also indicates that near but above the critical 
temperature T > T c , the pion mass vanishes in the sym- 
metric phase in general with m 2 ~ YP 11 ^ as S ~ m jjj, 
owing to wave-function renormalization |]l0f . Some of our 
previous arguments may apply provided that due care is 
paid to the pion dispersion relation with Q 2 — ► Q 2 ~ v |Q. 
Similar behavior has been noted at the mobility edge of 
a metal-insulator transition in metallic grains 
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How does the present chiral phase transition compare 
to the metal-insulator (MI) transition in QCD in vacuum 
we discussed recently in d=4 Qj? For comparison, we 
note that in the latter and for long times the quark return 
probability scales as |fl|,[12| 
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with a multifractal exponent jlj 
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7 = 2 X = 4- 
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In QCD j3 = 2 and the critical conductance er* is equal 
to the microscopic conductance ai — 2F 2 l 2 /ir for one 
mean- free path I, that is [fyjlq] 



(23) 
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where F w 93 MeV is the pion decay constant, E w (250) 
(MeV) 3 is the vacuum quark condensate and toq « 300 
MeV a typical constituent quark mass. Hence rj/4 w 
0.057, which yields a larger critical exponent than in (|is|). 
We note that our prediction for the critical exponent in 
QCD is in qualitative agreement with the numerical es- 
timates rj/4 w 0.08 — 0.16 obtained subsequently using 
the instanton liquid model |15| . Remarkably, the small- 
ness of 77/4 makes the result (pi]) close to a 2-dimensional 
diffusion process. 

In the MI transition, (^3|) follows from the ballistic re- 
gion and is related to a careful consideration of the null 
spectral sum rule pj JLTI , ^ . At the critical point this is 
necessary because of the multifractal character of the cor- 
responding wavefunctions at the mobility edge p2[ . In 
the metallic regime, the quark wavefunctions are spread 
through the metal, while in the insulator regime they are 
localized. At the edge, the wavefunctions are 'filamen- 
tary' with self-similar structure |rij| . In a second order 
transition we expect the quark wavefunctions to be per- 
sistently 'metallic'. 



V. PERCOLATION TRANSITION 

Multifractal wavefunctions have been deemed as 'quan- 
tum' percolation by Aharony Jig] . Could a 'classical' 
percolation transition take place at finite temperature in 
QCD? Recently, some arguments have been put forward 
by Satz ]l9|| . In the present framework this can be ad- 
dressed by noting that ([ll]) valid at Tc by standard scal- 
ing arguments |4|] does not incorporate the fact that for 
temperatures < T < Tc the Goldstone modes disperse 
asymmetrically in matter |2^j2l[]. Indeed, in these range 
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of temperatures and for space-like momenta we have in- 
stead of (H), 



2T,n 



DiQl 



D S Q 2 + 2m-i\ 



(24) 



with T,x = | (<ZS)t| the temperature dependent quark con- 
densate. Here, the 'temporal' D4 and 'spatial' D$ diffu- 
sion constants are related to the 'temporal' F4 and 'spa- 
tial' Fs weak decay constants of the pion ^0|. In the 
space-like regime they are both real, and we have gener- 
ically 



D 4 = 2F|/S T , 
D s = 2F i F s /X T . 



(25) 



If we were to denote by m\ the squared pion 'time-like' 
mass, and by m s the squared pion 'space-like' mass, then 
D4 = 2m /ml and Dg — 2m/m 2 s where m$ > m 4 by 
causality |20(| . Hence D4 > Ds- 

In terms of ( pi] ) the quark return probability in the 
range < T < T c reads 
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where we have used the facts that £t = t^Pt, At = 
1/prV and N — E/ At- At finite temperature both the 
density of states pr and the level spacing At change in 
a 4-volume V = L 3 /T. The momenta in ( p6f ) are Q4 = 
ni2irT and Q — n2n / L. ( ^6| ) implies the existence of two 
Thouless energies E4 — D4T 2 and Es — D$/L 2 , hence 
two ergodic times T4. crg = I/E4 and rg^rg = 1/-E-S, with 
T 4,crg < T s,crg- The corresponding Ohmic conductances 
are 174 = E^/At and 175 = Eg/ Ax- These scales allow 
a simple organization of the disorder in an asymmetric 
(finite temperature) Euclidean box. In particular, ( p6| ) 
becomes universal for mV <C 1 and times t ^> Tg )erg > 
T4 j01g . Most of the vacuum arguments presented in 0] 
can be extended to the present finite temperature phase. 

If we were to denote the 'temporal' £14 and 'spatial' gs 
quark conductivities, then by the Kubo formulae: g — 
DpT, we have g 4 > gs- For temperatures < T < Tc the 
quark conductivity in the 'spatial' directions is weaker 
than the quark conductivity in the 'temporal' direction. 
This is easily understood by the fact that 'spatially' the 
system 'screens' all charges including the singlet ones. 
Indeed, the pion screening length at high temperature 
asymptotes m| ~ 2nT by standard arguments ||, so 
that the spatial conductivity Ds ~ 2m/(27rT) 2 is para- 
metrically small in the high temperature phase. 

The asymmetry in the conduction properties may 
cause an 'asymmetric' percolation from d=4 to d=l as 
Ds becomes much smaller than D4, a situation reminis- 
cent of finite density ]22|j. From (123), it follows that 
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in the limit D4 3> Ds- We note with Aharony that in 
a 'classical' percolation transition in cither d=3 or d=4, 
the quark return probability on an infinite critical cluster 
is expected to scale as 
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where D — 2.5,3.2 is the fractal dimension and 9 = 
1.77,2.70 in d=3 and d=4 respectively. Hence, D/(2 + 
9) = 0.66, 0.68 for d=3,4 respectively, which are similar 
to 0.5 in (27). Both exponents are markedly different 
from the MI behavior (|2l|). 

In a linear sigma model the diffusion constants D4 and 
Ds can be estimated in a weak-coupling (1/n) expansion 
p0| . The results in our case and to 1-loop are 
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where m a ~ 500 MeV is a typical sigma mass. To the 
same-order we have St ~ S (1 — T 2 /8F 2 ). Similar calcu- 
lations can be performed in the context of the non- linear 
sigma model as well using instead 2-loo ps. A fine tuning 
of the running logarithms in this case pl| should bring 
both estimations into agreement. From (|2^) it follows 
that an order of magnitude drop in the relative conduc- 
tivities gs/ 04, ~ 0.1 or equivalently 
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(30) 



takes place for T ~ 180 MeV, which is were we expect 
the percolation transition to possibly set in. At this tem- 
perature only 50% of the chiral condensate or the density 
of quark states at zero virtuality has depleted by the 1- 
loop argument. It would be interesting to investigate the 
effects of strangeness on the present arguments. 



VI. SEMI-CLASSICAL APPROXIMATION 

What does the result (|ll]) mean for the light quark 
spectrum in QCD? A qualitative answer to this question 
can be obtained using semi-classical arguments p3| , p4j . 
Since in our case the density of states p(X) ~ A 1 /* 5 changes 
in the energy band of interest, the standard semi-classical 
arguments have to be modified. Indeed, a rerun of the 
standard arguments yields for the quark return 

probability 



P® = ¥ r ( J E\A j (x)\ 2 6{t-T j ) 
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and the spectral form-factor 



A,- (x) 
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(32) 
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where the sum is over closed classical paths of reduced ac- 
tion Sj(A), virtuality A, and period Tj — dSj/dA. Aj(x) 
is the Gaussian contribution around the classical paths 
to the propagator S(x,x;m — i(A ± A/2)). The aver- 
aging in (p2) is over the gauge-configurations (static or 
t-independent disorder). 

The contributions to ||l]) result from closed classical 
paths labeled by j, starting at x and returning back to x 
with probability Pj = \Aj(x)\ 2 . The contributions to j3^ ) 
differ from those of ( pl| ) by the overlap integration over 
the '4- volume' of a given classical path |2^,|24|]. For a 
virtuality Xj ~ 1/7} about the mean A ~ 0, the '4- 
volume' is of order Tj/p(l/Tj), as p(A) counts the num- 
ber of states per '4-volume' per virtuality. Inserting this 
result into (|32|) and using (|3l]), we conclude that 



A2|+|l+l/5 

*(*)«K-^p(t) 



(33) 



The dimensionless coefficient K is not fixed by these es- 
timates. For QCD = 2 §]. © in combination with 
(H) yields K(t) ~ |i| for large times (m = 0). This re- 
sult at T = T c is reminiscent of the result at T = and 
suggests once more that it is universal. Hence, the corre- 
sponding spectral rigidity is still logarithmic but with a 
coefficient that is not fixed by the present semi-classical 
arguments. It can be fixed using our recent arguments 
on critical scaling in 0-dimension [pl|2£ 



VII. CONCLUSIONS 

We have shown that the long time behavior of the 
quark return probability in a QCD phase transition car- 
ries a quantitative information on the character of a 
phase change, that could be used to discriminate be- 
tween a second order, a metal-insulator or a percolation 
transition. For times comparable to the Heisenberg time 
tji ~ l/A* and in the narrowing ergodic regime m <C A*, 
we have found that the quark return probability at the 
critical point scales as lil -1 / 5 for a second-order transi- 
tion. In contrast, in a metal- insulator transition it is of 
order |*| -0 - 943 while for an asymmetric percolation tran- 
sition it is of order |*| -1 ^ 2 from d=4 to d=l. Our results 
can be tested by numerically assessing (^) using current 
lattice QCD algorithms. 
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